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ABSTRACT 

Mixtures  of  inert  gcses  con  be  used  to  improve  performance  in  closed 
gos  turbine  cycles.  In  the  present  work,  neat  transfer  and  wall  friction  para- 
meters have  been  obtained  numerically  to  demonstrate  the  effects  of  mixture 
composition  and  gcs  property  variation  for  heating  or  coaling  in  regenerative 
heat  exchcngers  of  such  cycles;  the  situation  is  modelled  by  laminar  flow 
through  short  ducts  with  constant  wall  heat  flux.  For  design  predictions  ac- 
counting for  the  effect  of  property  variation,  it  is  found  that  the  property 
ratio  method  is  better  than  the  film  temperature  method  for  heat  transfer 
while  the  latter  method  is  preferable  for  apparent  wall  friction  - with  the 
proviso  that  the  present  definitions  of  the  .non-dimensional  parameters  be 
employed . 
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XI 


I 


nomenclature 


o,b,c 


Dh 

G 

9c 

h 

k 

L 


exponents  for  temperature-dependence  of  properties, 
equation  (6) 

specific  heat  at  constcnt  pressure 

hydraulic  diameter  (twice  plate  spacing) 

average  mass  flux 

dimensional  constant 

heat  transfer  coefficient 

thermal  conductivity 

length 

pressure 

exponents  for  property  ratio  method,  equations  (llj/and  (15) 

absolute  temperature  / 

* 

axial  velocity 
bulk  velocity 
axial  coordinate 


Greek  symbols 

force  constants  in  lennard-Jones  (6-12)  potential 
JJ  absolute  viscosity 

p density 


Dimensionless  quantities 

f'  fap  apparent  friction  factor  based  on  one-dimensional  apparent 

wall  shear  stress,  e .g.,  equation  (13);  f /lengthwise  mean 
apparent  friction  factor  a 

enthalpy,  (H  - Ho)/(CpoTo) 

L length  for  velocity  boundary  layer,  4L/(D^Re) 

l length  for  thermal  boundary  layer,  4l/(D^RePr) 

Nu  Nusselt  number,  hDjy/k 


IV 


p 

Pr 

Q+ 

Re 

+ 

v 

+ 

X 

m 

x 

y 

Subscripts 

a 

b 

ep 

e 

f 

fd 

m 

ref 

w • 
x 

o 


2 

pressure,  2g=  PQp/G 

Prandtl  number,  c 

wall  heat  flux,  q"^D^/(l<aT  ) 

Reynolds  number,  GD.  Re  based  on  axial  coordinate 
transverse  velocity,  HAJvRe/V^ 

axial  coordinate  for  velocity  boundary  layer,  4 x/(D,  Re) 
axial  coordinate  for  thermal  boundary  Icyer,  4x/(D^R ePr) 
transverse  coordinate,  y/D^ 

overall  average  or  mean  value  (lengthwise) 
properties  evaluated  at  bulk  temperature 
based  on  constant  property  idealizations 
in  entry  region^ 

based  an  film  temperature,  (T^  + J^)/2 
fully  established  or  asymptotic  value 
lengthwise  mean  value 
reference 

properties  evaluated  at  wall  temperature 
properties  evaluated  at  local  temperature 
inlet  conditions 


The  circumflex  ( ) represents  non-dimensionalization  with  respect  to  the 

value  of  the  quantity  at  the  inlet,  e.g.,  p - p/pQ. 
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The  closed  3rayton  cycle,  or  gcs  turbine  cycle,  hcs  been  suggested  cs 
an  efficient,  comoacf,  versatile  system  for  power  plant  end  prcoulsion  appli- 
cations [l,  2j  . With  thermodyncmic  cycle  studies,  3cmmert  and  Klein  [3] 
showed  that  considerable  savings  in  the  tofcl  casts  of  a gas  turbine  cycle  can 
be  achieved  by  mixing  helium  with  a heavier  ges;  the  increase  in  density  re- 
duces the  size  of  the  turbomachines  while  the  reduction  in  thermal  conductivity 
increases  the  size  of  the  heat  exchangers  so  ‘hot  an  optimum  occurs  at  an  Inter - 
mediate  molecular  weight.  Approximate  calculations  indicate  that  when  the 
heavier  gas  is  ano‘her  noble  gas  further  improvements  in  heat  trensfer  perfor- 
mance are  possible,  compared  to  pure  gases  at  the  same  pressure, temperature 
and  molecular  weight  [4]  . However,  the  geseous  data  and  correlations  for 
heat  exchangers  have  been  obtained  with  pure  gases,  primarily  air,  with  negli- 
gible vcriction  of  the  transport  properties.  Whether  such  results  con  be  applied 
for  mixtures  of  inert  gases  with  temperature-dependent  properties  is  a basic 
question  which  the  present  paper  attacks. 

In  gas  turbine  cycles  the  regenerative  heat  exchangir  is  typically  con- 
structed of  parallel  plates, with  short  fins  attached  forming  additional  parallel 
surfaces.  Consideration  of  the  heat  transfer  performance  versus  pumping  power 
requirements  of  these  heat  exchangers  usually  results  in  design  for  operation 
in  the  laminar  or  transitional  flow  regime.  For  leminar  flow  heat  exchangers, 
the  sfreamwise  length  of  the  fins  is  shortened  in  order  to  teke  advantage  of 
the  increased  heat  transfer  coefficient  of  developing  boundary  layers  by  conti- 
nually reinitiating  the  boundcry  layer.  The  thermal  boundary  condition  is  an 
approximately  constant  wall  heat  flux.  As  a guide  to  the  effects  of  mixture 
composition  and  property  variation  in  such  geometries,  the  present  peper  - for 
a first  objective  - investigates  the  simultaneous  development  of  leminar  thermal 
and  velocity  boundary  layers  in  the  entry  region  of  parallel  plate  ducts. 

The  following  section  briefly  discusses  related  work  which  can  be  ex- 
tended to  provide  improved  guidance  to  the  designer  of  regenerative  hear  ex- 
changers for  mixtures  of  noble  gases.  Section  3 summarizes  pertinent  knowledge 
of  their  transport  properties  and  demonstrates  the  generalizations  possible  to  re- 
duce the  analytical  task.  Since  tffii  consequent  governing  equations  are  non- 
linear and  coupled,  they  are  solved  numerically  as  outlined  in  section  4.  The 
results  of  interest  in  design -lengthwise  mean  parameters  - are  presented  in 


lection  5 : first,  the  effect  of  composition  at  low  heating  rates  and,  then, 
the  effects  fhcf  the  temperature-dependence  of  the  transport  properties  cause 
on  the  heat  transfer  parameters  and  the  wall  friction  parameters,  separately. 
Finally,  the  major  conclusions  are  reiterated  In  the  lest  section. 
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2.  PREVIOUS  WORK 

In  preliminary  design  studies  for  closed  gcs  turbine  systems,  Vcnco  [4j 
estimated  relctive  hect  transfer  coefficients  and  pressure  drop  for  bincry  mixtures 
of  the  inert  gases  with  geometry  held  constant.  For  heat  trcnsfer  he  employed 
a constant  orooerty  version  of  the  Sieder-Tcte  relationship  suggested  by  Kem 

[5] 


1.86  Re  Pr 


1/3  Q.  1/3 


(o/i: 


= 2.95  (L*1 


This  relation  is  essentially  a thermal  entry  correlation  of  the  leveque  form  [6j 
which  strictly  cpplies  only  for  a linear  velocity  profile  as  in  the  wall  region 
of  a fully  established  flow.  With  long  tubes  the  Nusselt  number  should  become 
constant  rather  than  tending  to  zero  as  In  the  above  relationship.  For  pressure 
drop  calculations,  Vanco  used  the  friction  factor  for  fully  developed  flow. 

His  approach  is  still  prevalant  in  industrial  design  of  lamincr  flow  hect  ex- 
changers . 

As  shown  later,  one  of  the  main  effects  of  varying  mixture  composi- 
tion is  to  vary  the  Prcndtl  number.  If  the  velocity  profile  is  fully  established 
before  heating  commences,  the  function  Nu  0.*)  should  depend  an  geometry 
ond  thermal  boundcry  condition  and  be  independent  of  Pr  since  the  non- 
dimensional  variables  can  be  defined  so  that  Pr  does  not  appear  in  the  gover- 
ning energy  equation  or  its  ooundary  conditions  [7]  . Since  thermal  end  shear 
boundary  layers  grow  at  different  rates  when  the  Prandtl  number  is  not  unity, 
the  solutions  to  the  simultaneous  entry  problem  will  very  with  Prandtl  number. 
In  a numerical  analysis  Hwang  and  Fan  [8;  have  shown  this  variation  to  be 
significant  at  low  values  of  L*  for  constant  wall  heat  flux  and  constant  fluid 
properties.  For  the  compcrable  constant  wall  temperature  problem,  Schlijnder 
[9j  suggests  applying  the  Pohlhausen  solution  in  the  immediate  entry  as 

Num  e = 0.66-4  (L*f,/2  ?r’,/6  (2) 


in  continuous  functions  of  the  form 


nl  nu" 


for  design  comoutation . 
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Schade  and  McEligof  [lc]  developed  a numerical  solution  to  examine 
the  effect  of  air  property  variation  for  the  simultaneous  entry  problem  with 
strong  hecting  or  cooling  cpolied  to  two  parallel  plctes.  They  concluded  thct 
the  local  Nusselt  number  Increases  slightly  and  local  friction  factors  increase 
severely  with  heating.  To  engineers  accustomed  to  constcnf  property  onclyses, 
the  occurance  of  a Icrge  chcnge  in  friction  factor  and  pressure  drop  while  the  Nusselt 
number  is  changing  only  slightly  might  be  surprising.  Schcde  and  A<\coligot  showed  that 
such  comparisons  are  sensitive  to  the  choice  of  the  temperature  at  which  the  properties 
in  the  non-dimensional  parameters  are  evaluated  end  that  use  of  the  constant  oro- 
perty  idealization  can  lecd  to  either  dangerous  or  conservative  design,  deoen- 
dtng  on  the  application. 

For  the  heat  exchanger  applications,  designers  employ  mecn  parameters 
end  resort  to  empirical  methods  to  correct  for  fluid  property  variation  ['ll]  . 

The  two  most  common  schemes  are  the  reference  tempercture  method  and  the 
property  ratio  method  [l2j  . Based  on  experiment  end  approximate  cnalysi:  , 

Kays  recommends  exponents  for  the  latter  method  for  verious  geometries  end 
fluids.  However,  despite  the  ready  availability  of  appropriate  numerical 
methods  for  laminar  flows  for  over  a decade  [l3],  these  empirical  correlations 
have  not  been  refined  end  the  question  as  to  which  is  the  mere  accurate 
method  has  not  been  answered.  Accordingly,  a second  objective  of  the  pre- 
sent paper  is  to  exemine  this  question  and,  if  possible,  to  improve  the  expo- 
nents used  in  the  property  ratio  method  - for  the  pcrallel  plate  geometry 
and  constant  wall  heat  flux. 
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3.  TRANSPORT  PROPERTIES  OF  NOBLE  C-AS  MIXTURES 

While  reel  gas  properties  con  be  employed  In  fhe  numerical  ana- 
lysis in  fabular  or  “auction  form  directly,  it  is  advantageous  to  exploit  the 
similarities  between  different  mixtures  to  reduce  the  number  of  computations 
necessary  to  cover  the  range  of  conditions  of  interest.  Thus,  if  their  behaviour 
can  be  generalized,  fewer  variations  of  parameters  are  required  and  the  results 
are  more  concise  and  have  greater  usefulness.  Accordingly,  this  section  sum- 
marizes our  knowlegde  of  the  pertinent  transport  properties.  In  conjunction 
with  the  following  section  it  demonstrates  that,  as  a first  approximation,  the 
variation  in  mixture  composition  can  be  reflected  in  the  anolysis  in  terns  of  a 
single  variable  parameter,  the  inlet  Prandtl  number. 


The  Lennard-Jones  (6-12)  potential  can  be  employed  in  fhe  Chcpmcn- 
Enskog  kinetic  theory  to  predict  thermal  conductivity,  viscosity  and  Prandtl 
number  of  binary  mixtures  of  fhe  inert  geses  [l4].  There  has  been  considerable 
experimental  study  of  the  pure  gases,  particularly  helium  end  argon,  but  few 
data  on  their  mixtures  exist  to  check  fhe  predictions. 

With  force  constants,  G’/'rt.  and  O , as  suggested  by  Hirschfelder , Curtiss 
and  3ird  [14]  fhe  prediction  of  fhe  helium  viscosity  falls  about  eight  percent 
below  fhe  data  of  Dcwe  and  Smith  [lo]  end  Kaleikar  and  Kesfln  [id]  ct 
temperatures  around  900°C . Likewise,  fhe  predicted  thermcl  conductivity  is 
about  nine  percent  lower  than  fhe  measurements  of  Saxena  and  Scxenc  [^7j 
up  to  1100°C.  Similar  discrepancies  occur  for  xenon,  but  agreement  with 
argon  data  is  close.  In  general,  agreement  is  good  near  room  temperature  for 
these  gases. 

Thornton  [l8j  measured  viscosity  of  fhe  helium/xenon  system  at  20°C . 

His  data  agree  with  fhe  predictions  to  within  a few  percent.  On  fhe  other 

hand,  fhe  thermal  conductivities  obtained  by  Meson  end  von  Ubish  [ 1 9]  at 
0 

520  C show  an  increasing  divergence  frem  fhe  predicted  curve  as  fhe  fraction 
of  helium  Is  Increased.  Using  force  constants  suggested  by  DiPiopo  and  <esfin 
,20j  ror  fhe  helium  component,  instead  of  these  at  Hirschfelder,  Curtiss  and 
Bird,  leads  to  essential  agreement  with  the  recommended  vclue  from  fhe  Thermo- 
physical Properties  Resecrch  Center  C 2 1 j . 


Thermal  conductivity  and  viscosity  ere  presented  egeinsf  logcrithmic 
coordinates  in  Figure  i For  helium,  xenon,  ergon  end  some  of  their  bincry 
mixtures.  These  values  ere  based  on  the  Lenncrd-Jones  '6-12)  potential  with 
the  force  constants  C = 2.158,  3.292  and  4.055  A,  £/7<  = 36.20,  152.75 
end  229  °K  for  helium,  argon  end  xenon,  respectively.  The  xenon  values  ere 
rrom  Hirschfelder,  Curtiss  end  3ird  end  the  others  from  DiPipoc  end  <esrin. 

Viscosities  of  the  noble  geses  and  their  mixtures  differ  only  slightly 
with  molecular  weight  (composition) . The  variation  with  temperature  is  appro- 
ximately the  seme  for  all.  Using  an  idealized  temperature  dependence, 

( >J/ ^fe.)  = (T/Tref  )Q  (4a) 

one  Finds  the  exponent  "a"  ranging  From  0.7  to  0.8.  On  the  other  hend,  the 
values  of  the  thermal  conductivity  vary  over  an  order  of  magnitude  from  xenon 
to  helium.  Again  the  temperature  dependence  is  about  the  same  for  ecch  and 
con  be  idealized  cs 

( k / k ref  ) = ( ^ ^ re  f ) ^ (4b) 

For  the  mixtures  the  exponent  "b"  Is  typicclly  slightly  less  then  "a"  end  renges 
from  about  0.7  to  0.75.  As  o first  approximation,  "a"  end  "b1"  could  be  fcken 
as  equal  and  the  same  vclue  could  be  used  for  ecch  of  the  mixtures.  The  spe- 
cific heat  is  Independent  of  tempercture  but  varies  with  composition . 

Perhaps  the  main  surprise  is  the  Prandtl  number  variation  of  the  mixtures. 
Argon/helium  end  xenon/liel  ium  are  shown  in  Figure  2;  the  other  mixtures  yield 
curves  of  the  same  shape.  The  temperature  dependence  is  almost  negligible 
since  the  power  law  exponents  for  k and  differ  so  little.  The  Prandtl  num- 
ber of  the  pure  gases  is  2/3.  However,  each  binary  system  shows  a minimum  ct 
intermediate  concentrations  (moleculcr  weight);  for  xenon/helium  it  is  Pr  *0.22 
at  room  temperature  and  is  particularly  broed.  Other  Prcndtl  number  minime 
are  : krypton/helium,  0.32;  argon/hel  ium,  0.42;  krypfon/neon,  0.53;  xenon/ 
argon,  0.57.  Thus,  Prandtl  numbers  in  the  renge  0.5  to  2/3  ccn  be  obtained 
with  several  choices  of  bincry  mixture  and  concentration . 


properties  of  ooble  gas  mixtures,  Numf»er  iodical 
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4.  ANALYSIS 

The  details  of  the  numerical  analysis  ere  straight  forward  and  effecti- 
vely the  seme  as  other  methods  for  solution  of  coupled  parabolic  partial  dif- 
ferential equations  (e.g.,  reference  [22]  ) so  only  the  essentials  will  be  out- 
lined here.  The  number  of  parameters  which  must  be  varied  to  describe  the 
range  of  conditions  possible  is  determined  by  examination  of  the  non-dimen- 
sional equations  governing  the  problem. 

Under  the  steady,  internal  boundary  layer  approximations  - plus  the 
assumptions  that  (a)  mixture  concentration  remains  constant,  (b)  Mach  number 
<<  (c)  RePr  5 100,  and  (d)  natural  convection  is  negligible  - the  gover- 

ning equations  can  be  written  : 


Continuity : 


9pu 
d x * 


9 P v 
3 y 


= 0 


(5a) 


x -momentum  : 


P u 


Energy  : 


3 u 

dx* 


3 H 


dy 


.♦  3H 


J.  AL._L  3 

2 d x 4 d y 


P u dr-7  ♦ p = 

OX*  d y 


Integral 
continuity  : 


4 Pr.  3 y 

J P u d y « -i- 


3 T 
3 y 


(5b) 


(5c) 


(5d) 


A circumflex  (A)  represents  non-dimensional ization  with  respect  to  the  value 
of  the  quantity  at  the  entrance. 


The  gas  properties  may  be  idealized  as  : 


P.  T 1 


Tai  &=Tb; 


«p  * T d 


(6) 


Initial  conditions  are  : f (0,?)  * I,  p specified  and  C (0 ,y)  * 1,  i.e.. 


uniform  entry.  3oundary  conditions  are  (1)  the  nonslip  condition  for  velocities, 
(2)  constant  wall  heat  flux. 


ond  (3)  symmetry  of  the  flow  and  boundary  conditions  with  respect  to  the  centre 
plane. 

Examination  of  equations  (5),  (6)  and  (7)  shows  the  free  parameters  of 

the  mathematical  statement  are  Pr^,  Q^",  p^  end  the  property  exponents  a,b 

and  d.  For  the  present  paper  a and  b ere  taken  as  0.75,  d as  zero  and  o’ 

o 

is  set  sufficiently  high  that  the  Mach  number  is  small  in  the  range  of  interest. 
The  parameters  PrQ  (corresponding  to  the  mixture  molecular  weight)  and  Q+, 
the  heating  rate,  remcin  variable. 

The  problem  is  solved  numerically  with  progrem  BAND,  developed  by 
Greif  and  McEligot  [23]  for  flows  between  parallel  plates  with  thermal  radi- 
ative interaction  using  a bond  cbsorptance  model.  For  the  present  calculations 
the  capability  to  handle  thermal  radiation  was  suppressed,  but  property  varia- 
tion was  included  by  choosing  non-zero  values  of  the  exponents  a end  b.  The 
numerical  program  is  a finite  control  volume  analysis  using  implicit  algebrcic 
equations  to  represent  the  governing  equations;  these  equations  are  itera- 
ted ot  each  axial  step  to  treat  their  coupling  and  the  nonlinear  terms. 

Mesh  spacing  increases  in  both  the  transverse  and  axial  directions.  For 
the  results  reported  here,  81  transverse  nodes  were  employed  with  the  first 
usually  at  (y/D^)  = 0.001,  and  longitudinally  there  wens  ?0  steps  per 
decade  normally  Jtarting  at  xg  3 10  ^ . As  noted  by  Worsoe-Schmidt  and 
Leppert  [l3]  the  boundary  layer  approximations  are  .not  appropriate  for 
x*  < 10'J  so  no  results  are  reported  for  the  initial  two  axial  decodes.  Prior 
tests  have  shown  convergence  within  2 percent  for  Nu  ond  within  about  1 
percent  for  f with  this  grid. 


5. 


RESULTS 


Predictions  have  been  obtained  for  the  rcnges  0.2  6 Pr  ^ 2/3/ correspon- 
ding to  the  Prandtl  number  variation  from  mixtures  to  pure  gcses,  and 

•2  S Q * ICO.  For  O ^ 10,  p was  token  os  10^  giving  cn  inlet  Mach 

^ 4*  — 4 

number  of  0.035  while  at  higher  O , p^  = 10  for  ~ 0.011  . Pertinent  wall 

parameters  ere  listed  in  Appendix  A for  all  conditions  studied  in  this  investi- 
gation . 


With  heating  at  a constant  wall  heat  flux,  the  local  bulk  tempercture 
increases  continuously;  with  constant  specific  hect  this  increase  is  linear.  The 
viscosity  increases  with  temperature,  so  the  local  bulk  Reynolds  number 
<GV-“bx>  decreases  in  the  axial  direction  end  the  flow  would  be  expected  to 
remain  laminar.  Density  is  inversely  proportional  to  temperature  so  the  flow 
accelerates  as  it  is  heated,  also  normally  a stabilizing  influence.  This  conti- 
nuous acceleration  prohibits  the  occurence  of  invariant  velocity  and  temperature 
profiles.  The  wall  temperature  is  larger  then  the  bulk  tempercture  in  order  to 
transfer  energy  to  the  gas  so,  at  a given  cross  section,  the  viscosity  and  thermal 
conductivity  will  be  higher  near  the  wall  and  the  density  will  be  lower.  Con- 
sequently, parameters  defined  in  terms  of  wall  properties  will  have  different 
values  than  those  using  bulk  properties  in  their  definition.  While  the  increase 
in  viscosity  is  expected  to  increase  the  wall  shear  stress  and  decrecse  the 
velocity  near  the  wall,  thus  increasing  thermal  resistance,  the  increase  in  ther- 
mal conductivity  is  a factor  tending  to  reduce  thermal  resistance.  Likewise,  the 
expansion  due  to  reduced  density  counteracts  the  decrease  due  to  viscosity  to 
some  extent.  Near  the  entrance,  the  increasedvisccsity  and  decreased  density 
at  the  wall  augment  the  non-slip  condition  causing  transverse  flow  awey  from 
the  wall;  this  flow  also  carries  thermal  energy.  Further  downstream  the  trans- 
verse flow  decreases.  (These  various  effects  are  reversed  with  cooling).  While 
the  individual  effects  of  these  phenomena  can  be  forecast  in  some  cases,  their 
combined  effect  is  not  obvious  and  may  vary  with  geometry.  Since  density, 
viscosity  and  thermal  conductivity  variations  are  all  of  the  same  order  of  mag- 
nitude, none  dominates  in  such  a way  that  a single -parameter,  closed-form  ana- 
lysis would  appear  fecsible. 

For  the  present  study  the  presentation  of  results  emphasizes  mean  para- 
meters which  are  useful  to  design  engineers;  however,  loccl  Mussel t numbers 
and  friction  factors  ore  included  in  the  Tables  for  those  interested.  With  pro- 
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p«rty  variation,  parameters  can  tcke  a number  of  values  depending  on  how 
the  temperatures  which  are  used  to  evaluate  the  pertinent  properties  ere  de- 
fined. To  avoid  later  confusion  in  the  use  of  these  predictions  it  is  necesscr/ 
to  be  explicit  in  the  definitions. 

Local  bulk  temperature,  T^,  is  the  temperature  corresponding  to  the 
total  enthalpy  flow  at  a cross  section,  i.e.,  the  socalled  mixing  cup  tempera- 
ture. Average  bulk  temperature,  is  the  arithmetic  averege  of  the  local 

value  at  length  L and  the  inlet  value 

Tba  ~ ^bo  + TbL)/2 

While  an  integral  average  Is  normally  used  for  the  average  wall  temperature 
In  analyses  it  is  not  of  use  to  the  designer  who  lacks  the  detailed  knowledge 
of  Tw(x);  accordingly,  we  define  as  an  arithmetic  average 

Twa  = ^wo  + Twl)/2 

, which  becomes,  for  constant  well  heat  flux, 

Twa  " ^bo  * TwL>/2 

(Thu*,  is  lower  than  /T^(x)  dx/L).  Then,  in  agreement  with  normal- prac- 
tise, an  overage  film  temperature  is  chosen  as 

Tf=  - ^wo  + Tbo>/2 

Comparable  subscripting  is  used  for  indentifying  the  temperature  ct  which  pro- 
perties are  evaluated,  anc*  ,’m^ar^y  ? or  non-dimensioncl 

parameters,  Re^  = GD^/,U,  . Average  pressure  is  also  token  at  the  arithme- 

tic average  of  inlet  and  exit,  L.  Other  quantities  will  be  defined  later  as  used. 

Predictions  for  constcnt  fluid  properties 

For  small  heating  rates  or  low  temperature  differences  the  constcnt  pro- 
perties idealization  is  often  adequete  for  practical  cppl ications . For  example, 
one  sees  in  Figure  1 that  a 50°C  difference  in  temperature  causes  a chcnge  in 
thermal  conductivity  of  about  17%  at  room  temperature  and  4 % at  1CC0°K. 

In  Figure  3 are  plotted  the  mecn  Nusselt  numbers  versus  non-dimensional  length, 
L * 4L/(D^RePr),  as  obtained  by  setting  the  property  variation  exponents,  a 
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and  b,  at  zero  and  bolding  p - 1 . The  mean  Nusselt  number  is  defined  in 
terms  of  the  integrcl  cverage  heat  trcnsfer  coefficient  as 


h a • 0( 

k 


r 

kL  J 


h ( x ) dx 


From  Figure  3 one  sees  that  in  the  entry  Nua(L*)  increases  as  the 
Prandtl  number  is  reduced  by  employing  different  inert  gcs  mixtures.  This  in- 
crease is  approximately  17%  as  Pr  varies  From  2/3  to  0.2.  Since,  for  con- 
stant properties,  the  development  of  the  sheer  boundary  layer  is  a function 
of  x*  = 4x/(D^Re)  and  is  independen<.of  Prandtl  number  (i.e.,if  is  a solution 
of  equations  (5a, b and  d),  the  velocity  profile  is  more  fully  developed  for 
Pr  - 2/3  then  for  Pr  = 0.2  at  any  specified  value  of  L*.  Consequently,  the 
increase  in  Nuq  coincides  with  a higher  average  velocity  gradient,  (3u/3y)Q, 
between  0 and  L*;  with  increased  velocities  near  the  wall  an  increase  in  heat 
transfer  parameters  is  reasonable. 


Also  plotted  on  Figure  3 is  the  Sieder-Tate  correlation  (1)  which  has 
been  employed  in,  practise  for  heat  transfer  calculations  in  comparable  situa-  • 
tions.  In  addition  to  differing  by  as  much  cs  a factor  of  two  in  the  range  of 
interest,  this  calculation  fails  to  show  the  proper  trend  with  L*  except  appro- 
ximately in  the  limited  range  0.01  L**  0.1.  At  L*  *•  0.001  the  mean 

Nusselt  number  varies  as  rather  than  as  in  the  Sieder- 

Tate  correlation.  Failure  of  this  correlation  to  account  for  Prandtl  number 
dependence  has  been  mentioned  earlier. 

Based  on  an  integral/superposition  method  Kays  [l2J  suggests  that  the 
local  Nusselt  number  for  the  simultaneous  growth  of  leminar  external  boundary 
layers  can  be  predicted  by 

Nu  = 0.453  Pr1/3  Re  1//2  (8) 

X X 

for  the  thermal  boundary  condition  of  a constant  wall  heat  flux.  For  the  im- 
mediate entry  where  T^x  es  Tq  and  *»  Vq,  this  relation  can  be  transformed 
to  a mean  Nusselt  number 

Num  # * 0.906  Pr’/3  ReQ  J/2Dh1/2A 1/2  (9) 

In  this  case,  equetion  (3),  the  form  suggested  by  Schltjnder  ; 9],  would  become 
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1 /2 

Nua  = [8.2352  + 1.3122/(Pr1/3L*)]  (10) 

This  equation  is  shown  on  Figure  3 for  c pure  gas,  Pr  = 2/3, end  con  be  seen 
to  agree  with  *Le  numerical  prediction  within  about  ten  percent.  For  the  mix- 
tures agree~ier1’  is  as  good  or  better;  the  comparisons  are  shown  individually 
later. 


It  should  be  reemphasized  that  while  the  increcse  in  Nun.*)  is  of  the 
order  of  15  percent,  when  replacing  a pure  gas  by  a xenon-helium  mixture, 
the  gain  in  heat  transfer  coefficient  is  much  greater.  For  example,  taking 
equation  (9)  as  an  approximation  for  short  fins  in  conjunction  with  Figures  1 
and  2,  one  can  see  that  for  a given  geometry  and  Reynolds  number  the  effect 
of  replacing  pure  argon  by  a helium/xenon  mixture  is  to  increcse  the  heat 
transfer  coefficient  about  2.4  times. 

Predictions  of  heat  transfer  with  property  varioticn 

When  fluid  properties  very  significantly  due  to  high  heating  rates  and 
related  large  temperature  variations  in  the  flow  field,  the  numerical  values  of 
non-dimensional  parameters  such  cs  Nu  and  Re  depend  on  the  temperature  at 
which  their  properties  are  evaluated.  Two  methods  of  accounting  for  "the  pro- 
perty variation  are  common  in  practise  : the  film  temperature  approach  and  the 
property-ratio  approach  [12]  . In  the  fil  m temperature  approach  it  is  assumed 
that  using  properties  evaluated  at  T^  will  allow  direct  use  of  predictions  ob- 
tained under  the  constant  property  idealization.  In  the  property  ratio  approach 
the  bulk  temperature  is  used  for  the  properties  end  the  effect  of  heating  is 
represented  as 

Nu/NUcp  = tryyp  (1!) 

for  gases.  The  effectiveness  of  these  two  methods  for  heat  transfer  predictions 
is  examined  in  the  present  section;  wall  friction  predictions  are  considered  in 
the  following  section . 

Figure  4 demonstrates  the  epparent  effect  of  heating  rate  on  the  mecn 
Nusselt  number  when  the  properties  ere  token  at  T.  / i.e.,  the  property  ratio 
approach.  Thus,  l*^  is  defined  as  ^/(D^Re^Fr)  * 4L  kfcg/fCpGD^2)  • 
effect  of  heating  the  gas  in  0 slight  Increase  is  Nu^  compared  to  the  predic- 
tion of  Nu  for  constant  properties  (Q4"-*0)  at  the  same  value  of  l^.  For  a 
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given  value  of  O , Nuj_a  's  increased  more  for  the  pure  gases  than  for  the 
mixtures,  but  this  result  is  partially  a consequence  of  the  lower  wall-fo-bulk 
terpperature  reached  by  the  gcs  with  the  lower  Prandtl  number.  However,  for 
Q = 100,  which  leads  to  (T  /T,  ) around  2.5  to  2.7,  the  increase  is 
only  of  the  order  of  five  percent.  Thus,  the  exponent  p in  ecuction  (11)  would 
be  about  0.05,  which  ccn  probably  be  considered  negligible  for  most  prccfical 
cases.  On  the  other  hcnd.  For  Q =100  and  Pr  = 2/3  the  local  Nusself  number 
Nu^x  reaches  a peck  increase  over  the  constant  property  prediction  of  about 
twenty-five  percent  at  a location  where  » 3 so  the  meon  Nusself 

number  is  affected  even  less  than  the  local  value. 


To  examine  whether  the  film  temperature  aocroacn  is  cn  improvement, 
the  results  are  plotted  in  Figure  5 with  prooerties  evaluated  at  the  avercc° 
film  temperature,  Nu^  - end  = ‘Hkp^c^GD^ . With  this  choice 

the  mean  Nusselt  numbers  decrease  os  the  heating  rate  Is  raised.  This  effect 
is  of  approximately  the  seme  magnitude  for  the  pure  geses  as  for  the  mixtures 
wi  th  their  I ower  Prandtl  numbers.  The  most  significant  observation  is  fhet  the 
reduction  is  about  three  times  as  great  cs  the  change  when  using  bulk  proper- 
ties (Figure  4)  so  there  is  no  advantage  in  accuracy  with  the  film  temperature 
approach . 


in  me 


When  heafmg  fhe  gas  and  bofh  Increase,  more  rapidly  Ir 
entrance  region  and  then  both  at  approximately  the  seme  rate  at  larger  distan- 
ces. Consequently,  first  increases  and  then  approaches  unity  axiclly. 

Likewise,  Nu^O-*)  first  greater  than  Nuc^(L*)  and  then  approaches  Nuc^(l*) 
downstream.  In  this  situation  correlations  of  the  form  of  equation  (11)  make 
sense.  In  contrast,  when  cooling  fhe  gcs  - as  on  fhe  apposite  side  of  a re- 
generative heat  exchanger  - both  temperctures  decrease  continuously  such  that, 
in  the  limit,  Ty/rfcj  approaches  zero  or  yT w approaches  infinity.  The  rela- 
tive property  variation  ccross  the  flow  becomes  greater  downstreem  instead  of 
less  as  in  fhe  heating  situation . This  limit  is  not  likely  to  be  of  concern  in  a 
gas  turbine  cycle  since  fhe  lowest  femperefure  is  set  by  fhe  inlet  temperature 
of  fhe  side  to  be  heated,  but  might  be  a difficulty  in  cryogenic  acpl  leaf  ions . 

In  the  cooling  range  covered  in  fhe  scope  of  fhe  present  peper,  Nu,3a(L£a  ) 
remained  sufficiently  close  to  Nu^^fL*)  to  neglect  fhe  difference. 


If  one  calls  ecuction  (10)  fhe  Schl(jnder-<cys  relation  end  plots  It  cs 
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a dashed  line  on  each  of  the  subfigures  of  Figure  4,  if  may  be  seen  fhaf  agree- 
ment with  the  constcnt  property  result  is  goad  for  ecch  cose  but  successively 
better  as  Pr  is  reduced.  Since  Nu^  (L^)  is  predicted  so  closely  by  Nu_  'l*) 
it  is  worthwhile  to  improve  the  equation  so  if  mcv  be  used  by  the  designer 
over  the  range  of  interest.  By  adjusting  the  Prandfl  number  dependence  of 
Num  e for  better  prediction  at  1*  = 0.001,  one  may  obfcin 

Nuq  = [ 3 . 2352  + 1 .9312/(Pr0,25V  (12) 

This  equation  still  is  of  the  order  of  five  percent  lower  than  the  numerical 
prediction  near  L = 0.01.  For  moderate  necting  rcfes  if  would  be  within 
about  ten  percent  of  the  numerical  results  and  would  be  low;  in  heat  exchan- 
ger design  this  would  lecd  to  units  slightly  longer  than  necesscr/.  Alternatively, 
the  constants  in  equation  (12)  could  be  optimized  for  another  rcnge  at  the 
expense  of  the  accuracy  of  predictions  in  the  immediate  entry. 


Prediction  of  wall  friction  with  property  variation 

While  most  cnclyses  presently  available  for  developing  flows  present  fric- 
tion results  in  terms  of  the  wall  sheer  stress  evaluated  from  the  velocity  gradi- 
ent at  the  wall  (f$),  this  cpproach  is  not  of  use  to  the  designer  when  the  ve- 
locity profile  is  changing  substantially,  as  in  the  entry  or  when  a gas  is  heated. 
To  predict  the  required  pressure  drop  with  a one-dimensional  design  procedure, 
one  uses  the  "apparent"  friction  factor,  f , besed  on  the  wall  shear  determined 
by  treating  the  momentum  change  as  one-dimensional  . The  same  treatment  is 
often  employed  in  experiments  where  size  prohibits  velocity  profile  measure- 
ments. Both  methods  of  presentation  can  be  chosen  with  numerical  results;  con- 
sequently, Bcnkston  and  McEligot  t22j  were  able  to  demonstrate  (a)  the  nume- 
rical values  of  f and  f con  differ  substantially  and  (b)  discreoancies  ecrlier 
s ap  ' 

thought  to  exist  between  experiments  end  analyses  were  primarily  due  to  the 
differences  in  the  definition  used  for  the  friction  factors. 


As  with  the  heat  trensfer  results  we  concentrate  in  presenting  a mean 
copcrent  friction  factor, 


fQ  1 
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(13) 


(The  local  opperent  friction  factor,  f , appearing  in  the  Tables  Is  defined 
in  ffa#  craiogou*  dcrivotiv®  form  with  d/dx  replacing  (^  ) ).  When  P ^ 
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is  constant,  the  second  term  in  brackets  does  not  change  and  the  definition 
reduces  to  that  of  Shch  and  London  ^7j  . With  constant  fluid  properties  one 
solves  equct'or s (5c),  (5b)  and  ’5d)  only,  sc  the  result  is  indeaendent  of 
Prcndtl  numoer  and  con  be  written  cs  a single  function  f (L  ) which  ccproc- 

che>  fQ  • Re/ 24  as  L secomes  Icrge.  This  function  may  be  found  robulcted  in 
Aopendhx  A or  ccn  be  derived  from  earlier  loccl  results  (10,  24  ] . For  a con- 
tinuous approximation,  the  approach  of  Schlunaer  ccn  be  used  cs  in  equcficn 
(13)  to  give 

fa  • Rs/24  = \J  1 + 0.0788/L"'  (14) 

which  represents  the  numerical  results  well  in  the  immediate  entry  but  is  4 to 
5 percent  high  in  the  rcnge,  0.05  < L < 0.2. 

With  varying  transport  properties,  the  energy  equction  (5c)  is  coupled 
to  equations  (5a,  b end  d)  via  the  temperature-dependent  viscosity  end  den- 
sity, so  the  wall  friction  also  becomes  a function  of  the  Prcndtl  number  end 
the  heating  rate.  Again  the  question  arises  as  to  the  better  method  of  accoun- 
ting for  the  fluid  property  variations.  Predictions  of  friction  ere  not  cs  well 
behaved  cs-hecf  transfer  parameters.  In  contrast  to  the  heat  trensfer  results, 
direc*  use  of  the  overoge  bulk  properties  in  f • Re  end  L+  does  not  collcpse 
the  results  nicely  around  the  prediction  based  on  constant  properties;  the 

main  effect  is  to  spread  the  curves  towards  larger  l,  as  Q*  increases. 

ba 

The  effect  of  heating  rcte  on  apparent  well  friction  is  presented  in 
Figure  6 partially  in  terms  of  average  bulk  preoerties.  That  is,  P,  is  used 

SO 

for  the  coefficient  in  equation  (13)  end  Re,  is  defined  as  before -but  the 
non-dimensional  length  is  based  on  Inlet  properties,  i.e.,  lj".  With  this 
representation  heating  increases  f^Q  • Re^  considerably  more  than  Nu.  is 
raised  at  the  same  level  of  Q+  . At  lengths  greater  than  L+  = 0.1  the  curves 
with  heating  approach  the  constant  orcoerties  curve  only  slowly,  although 
T^^/T^  lJ  close  to  unity,  as  the  heated  entry  continues  to  affect  the  inte- 
grated results  far  downstream.  Close  inspection  of  the  trends  for  the  highest 
heating  rates  shows  thet  as  Pr  increases  c convergence  - crom  heefea  entry 
behavior  towards  agreement  with  constant  property  behcvlcr  - is  moved  further 
downstream.  This  effect  corresponds  to  the  difference  in  growth  of  the  thermal 
boundcry  layer  and  sheer  bounder/  layer  as  the  Prandf!  number  chcnges  : 

Nu(l*)  shows  only  a moderate  effect  or  Pr,  so  for  the  same  hearing  rate  khe 
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value  of  T^A^  is  clmost  the  seme  at  eque!  value  of  L*  rcther  then  l+,  thus 
T^/T^  approaches  unity  for  Pr  = 0.2  at  earlier  vclues  of  L+  then  for  Pr  = 2/3 
end  the  veriafion  of  properties  across  the  channel  is  less  for  Pr  = 0.2  at  the 
same  L . While  the  friction  predictions  for  heating  approcch  the  adicbcfic 
prediction  cs  L^a+  increcses,  those  for  coaling  diverge;  this  result  also  corres- 
ponds to  the  trend  of  property  variation  since  the  rctio  T.  /j  increcses 

□CT  wO 

downstreem  for  cooling  cs  described  earlier  in  the  section  on  heat  trensfer. 

As  with  the  heat  trensfer  results,  the  epparent  effect  of  property  vcric- 
tion  on  well  friction  is  sensitive  to  the  choice  of  reference  temperature.  With 
average  film  temperature  for  the  reference,  the  shape  of  the  resulting  curves 
differs  from, the  shape  with  bulk  temperature  as  reference.  In  Figure  7 the  pro- 
duct fpa  • R ®fa/2‘i  is  plotted  against  Lq+;  p p is  used  in  the  coefficient  in 
equation  (13)  and  Re.^  is  besed  There  is  no  advantage  in  comparison 

on  the  besis  of  LpQ  since  results  are  shifted  then  further  to  the  right  (with 
hecting)  so  thet  for  l+  > 0.01  the  difference  from  the  cdicbatic  prediction 
is  increased.  For  heating:  the  friction  parameters  ere  reduced  for  short  lengths; 
then  the  predictions  coverge  with  end  cross  the  ccnstcnf  properties  curve  and 

0 

remain  slightly  greater  at  larger  distances.  In  comparison  to  the  bulk  property 
predict  ions,  the  effects  for  short  and' long  ducts  are  approximately  the  same 
magnitude  with  strong  heating,  but  for  intermediate  lengths  and  for  G+  "3?  10 
a display  In  terms  of  film  properties  shows  significantly  less  variation.  In  the 
range  -2  < Q <2  there  is  no  significant  effect  of  heating  until  lQ+  approa- 
ches 0.1  with  film  properties  and  then  the  effect  is  only  of  the  order  of 
five  to  ten  percent.  As  is  the  case  for  bulk  properties,  the  convergence  to- 
wards the  adiabatic  prediction  is  at  successively  greater  distances  (lg+)  cs  the 
Prandtl  number  increcses,  but  for  the  seme  condition  it  is  several  times  eerlier 
with  film  properties . With  cooling  ; the  direction  of  the  trends  ere  reversed 
but  for  Q > -2  they  are  essentially  again  negligible  for  entry  problems. 

It  is  not  clear  from  Figure  5 end  7 which  approach  is  better  ; property 
rctio  or  film  temperefure.  The  property  ratio  approach  would  be  reoresented 

as 

<fbo  • R*0)/  <fo  • *a.X  * (Tw/rba)q  05) 

so  this  quotient  is  plotted  versus  temperature  ratio  in  Figure  3 to  examine  the 
suitability  of  a single  exponent.  A complicated  pattern  appears.  In  contrast 
to  the  expectation  of  Kays  and  London  [ 25;  , the  general  trend  Is  a substantial 


Bulk  properties 


proper  tier 


Examination  of  property  ratio  method  for  mean  wall  fricti 
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increase  with  temperature  ratio.  There  ere  slight  differences  with  Prcndtl  num- 
ber but  the  trends  are  mostly  the  same.  An  exponent  q of  the  order  of  unity 
would  overpredict  the  friction  factor  at  the  higher  ferroercture  reties  and  under- 
predict it  at  lower  values.  For  cooling,  q = 1 is  velid  within  a few  percent. 
For  moderate  heating,  the  necessary  value  of  q (i.e.,  the  slooe  of  a line 
from  the  origin  on  this  logcrithmic  plot)  varies  with  length  l*:  it  is  appro- 
ximately constant  as  the  temperature  ratio  increases  with  length  then  increases 
gradually  as  the  ratio  drops  for  successively  longer  ducts.  The  latter  effect 
is  a consequence  of  the  slow  covergence  of  f^Q  • Re^  to  the  adiabatic  curve 
for  long  ducts  as  discussed  earlier.  If  is  seen  that  a function  q (Lo+,  Q+,  Pr) 
would  be  necesscry  to  describe  the  derailed  behavior.  For  Q+  = 2 and 
*-0  <0*^/  exponent  q * 1 .5  would  reduce  the  difference  from  the  constant 


properties  curve  from  13  percent  to  a 7 percent  discrepancy.  With  Q4-  = 10, 

^ “ 1 *2  is  a better  approximation,  but  the  discrepancy  would  still  recch 
twenty  percent.  These  comments  and  comparison  of  Figure  6 and  S suggest 
that  the  two  methods  hove  approximately  the  same  overall  cccuracy  for  Q+  < 2 
with  a slight  advantage  to  the  Rim  properties  approach  for  short  ducts.  For  mo- 
derate heating  - to  Q = 10  - the  film  property  method  is  clearly  superior, 
while  at  higher  heating  rates  both  methods  show  regions  where  the  simple 
correlations  would  mislead  the  designer  substantially. 


If  is  perhaps  inconvenient  for  the  designer  to  have  one  method  perform 
better  for  heat  transfer  while  the  other  is  preferable  for  wall  friction,  but  the 
difficulty  should  be  negligible  provided  the  present  definitions  of  the  para- 
meters are  used.  Once  the  heat  transfer  problem  is  solved  for  the  wall  tempe- 
rature using  average  bul'x  properties,  the  average  film  temperature  can  be 
calculated  from  the  results  ana'  can  then  be  employed  to  predict  the  wall 
friction  behavior. 


Analytical  correlations  such  as  equations  (12)  and  (U)  are  useful  for 
parameter  studies  of  systems  and  for  initial  sizing  of  components  when  hundreds 
to  thousands  of  individual  configurations  may  be  calculated.  When  greater  ac- 
curacy is  needed  m Rnal  design  decisions  - or  if  variable  wall  hect  flux  should 
be  treated  - the  numerical  analysis  can  be  employed  directly.  With  the  direct 
application  of  the  program,  the  question  of  definitions  of  the  non-dimensional 
parameters  is  avoided;  the  engineer  con  choose  definitions  to  suit  his  own 
convenience,  including  direct  presentation  in  temperatures,  pressure  and  lengths 
in  units  of  his  choice. 
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6.  CONCLUSIONS 


For  Heat  transfer  to  mixtures  of  inert  gases  in  short  ducts  formed  of 
pcrallel  plates,  the  following  conclusions  concerning  the  behcvior  of  the  mean 
parameters  are  wcrranted.  In  terms  of  approoriate  non-dimensional  variables 
and  parameters,  the  effects  of  varying  mixture  concentration  can  be  represen- 
ted by  variation  of  the  inlet  Prcndtl  number.  These  conclusions  are  based  on 
the  specific  definition  of  parameters  chosen  in  the  present  work;  with  strong 
heating  rates  the  use  of  alternate  definitions  in  the  correlations  can  cause  sub- 
stantially different  predictions  of  the  heat  transfer  coefficient  and  friction 
factor  and,  consequently,  of  the  wall  temperature  end  pressure  drop. 


a) 


b) 


c) 


Under  the  constcnt  property  ideal  Izcfion,  heat  transfer  and  epoarent 
wall  friction  parameters  ccn  be  approximated  by 


Nu  * [a . 2352  + 1 .9312/(Pr°-254L')] 


1/2 


and 


1/2 


(fa  • Re/24)  = [ 1 + 0 .0783/l+  ] 
to  within  ten  percent  for  L*  and  L+,  respectively,  greater  than  0.C01. 


For  heat  transfer  in  the  range  -2  < Q+ < .100  the  bulk  properties/pro- 
perty ratio  method  of  accounting  for  the  effects  of  gas  property  vcric- 
tion  - with  exponent  p 5 0.005  - provides  better  predictions  than  the 
film  temperature  method. 


For  wall  friction  in  the  range  -2  <Q'  ^10  the  film  temperature  method 
is  more  accurate  overall  than  the  property  ratio  approach. 
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APPENDIX  A 
TABULATED  RESULTS 


This  appendix  presents  reproductions  of  the  summary  tables  of  the  com- 
puted results  that  were  provided  as  output  from  the  UNIVAC  1108  computer 
at  the  Rechenzentrum,  Universifdt  Karlsruhe.  All  cases  for  which  predictions 
were  calculated  are  included  in  these  tables.  Space  limitations  prohibit 
listing  all  the  predicted  parameters  of  interest;  however,  the  set  shown 
should  provide  the  basic  infomnation  from  which  others  ccn  be  ca'culated 
via  their  definitions  and  the  idealized  property  relationships.  For  example. 


LbaS  ^Re^Pr 


& l R*0 

DhRe0pr  Reba 


*>*  ( Tba  \° 

DhPfo  V T0  ) 


The  cases  calculated  are: 


Pr 


Q+ 


0.7  (air) 
0.666  ... 
0.4 
0.2 


50,  0+ 

100,  50,  25,  10,  2,  0,  -2 

100,  10,  Ot 

100,  50*  25*,  10,  2,  0, 

-0.2,  -2 


■‘■The  case  of  Q - 0 represents  heating  rates  sufficiently  small  that 
the  properties  can  be  considered  constant.  In  order  to' calculate 
these  cases  non- zero  heat  flux  was  chosen  and  the  property  exponents 
were  chosen  for  constant  properties;  these  are  identified'  in  the 
following  tables  by  the  head:  VTS  • .000,  CON  - .000,  CP  - .000. 

•Inlet  pressure  relatively  low,  leading  to  excessive  Mach  number  at 
last  few  stations. 
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Of  necessity  W'ne  of  the  fable  Heading*  are  brief  so  the  Following 
definitions  and/or  comments  a*e  appropriate  : 


Program 

Symbol 

Usual 

Nomenclature 

Definition 

Comment 

PR,  0 

Pr 

0 

(%fA>0 

Inlet  Prandfl 
number 

QPP-DH/KO-TO 

Q+ 

GAMMA 

> 

u 

\ 

u 

Specific  Hea* 
ratio 

VIS 

a 

viscosity  exponent  in 
equation  (6) 

CON 

b 

thermal  conductivity 
exponent  in  equation  (6) 

CP 

d 

specific  heat  exponent 
in  equation  (6) 

4X/DPE0 

• 

X 

o 

4x/(DhReQPr) 

TB/Tl 

Tb/ro  • ‘ 

local  bulk  temperature 

TW/TB 

Twx/Tbx 

local  temperature  ratio 

NUBX 

Nv 

hDf/Sx 

local  bulk  Nusselt  number 

FRB/24 

fbx’R*bx^4 

local  apparent  friction 
product,  based  on  bulk 
properties,  see  pg.14 

WA/BA 

Tw^ba 

average  femoenafure  ratio, 
see  pg . 10 

NUBA 

Z 

c 

8* 

W'ba 

avercge  bulk  Nusselt 
number,  see  pg.  11 

FRBA/24 

fba,R#b</24 

mean  apparent  friction 
product,  bcsed  on  bulk 

PO-PI 

K 


P "P 

rO  X 


properties,  eq.(13)  and 
pg.  15 


2gcPQfp3-p)/G' 


( Po  " Pu)  ” ^Pfca 

PbL 

2Ux, 


pressure  droo,  note  og.9 

A 

c ? 6a 


2-a  W3  Lv-1  J 

comparable  to  K(x)  c f 
Shah  and  London  £7]  - 
but  error  in  output  in 
earlier  constant  property 
results. 


i_A8X.OA*  (iAS.FLO'J  BET^EE'l  ^i°ALLEi  PLATES 
tlaEJM  tC  C3  4Sta.’.T  AALL  r»FAT  FLaA 


P9.1*  »73?»  0>.0»  Aa  \ • *<00 

®*<;,PE8t7  EXPONENTS,  V!a»  .a73,C0M»  .53r.CF»  .J?5 
8X/.JPE0  TO/  7 I r*/T^  H >J  R A r ^ 8/  2**  «a/RA  -iU^A  F93A/2H  PT-PI 


I 

• os 

2 

2 6 

36 

3a 

9 

• 235 

:-T.: 

. j 4 

2 

37  1 

32 

5 i 

7 

.7  75 
O e 3 

t . 

1 

• LL 

2 

2 

‘tT2 

ZlA 

2 9 
2S 

1 3 

a ^ a 3 

.112 

I • 

ijt 

I 

l 

* i s 
. ; •* 

2 

2 

a n O 

73? 

2 3 
2 1 

i ^ 

r3 

z 

5 

.157 

I * 

A • 

i 

i 

• - 2 

.27 

2 

2 

7 / 
7i3 

1 9 

1 3 

n * 
\ 4 

H 

9 

.39  4 

1 * 

2 A 

l 

1 

. J 7 
. 39 

2 

2 

7** 

t 6 

t S 

7 

9 

3 

..  u H ^ 
.789 

C -a  C 

4 « 
2. 

. 4 3 

2 

?. 

7:;3 

4*>9 

i* 

12 

>r» 

g7 

3 

3 

A O 2^ 
• 2 4 0 

Z . 

1 • 

l 

z 

» ^‘t 

• 1 3 

2 

2 

•4  J | 

739 

1 ! 
1 2 

V 

2 

7 

A 7 3 J 

,699 

I • 
1 • 

2 

2 

•7_« 

2 

1 

n2« 

A 3 u 

9 

S l 
87 

? 

2. 

a 3 *4  6 

• 1 l!L 

l • 

! r* 

3 

3 

• 1 1 
.ss 

1 

1 

7 1 2 
5 7 2 

9 

9 

*3  f 
8 7 

1 

1 

• a 7 4 
.672 

1 * 
l* 

-♦ 

C 

* ?•» 

.65 

_ A _ 

t 

1 

34  1 

3 6 7 
. 7 

9 

8 

7 4 
6 1 
c r. 

t 

a 5 7 J 
• * ) 2 

l ■ 

1 • 

9 

°3 
. 6 3 

t 

1 

Il-ii. 
1 25 

a u 

9 

rt 

“ -4 

3c 

A 

i 

.339 

. 1 99 

i ♦ 
1 • 

1 3 
1 7 
2 1 

. h J 

. 8S 
.22 

1 

1 

1 

r9» 

"35 

029 

8 

8 

8 

3-t 
2 1 
2? 

• 

1 

! 

.12' 

. )8» 

• O7  2 

t • 

1 . 

1 . 

645 
7 -i  7 
7?7 
3h  I 

7 3 1 
736 

979 
VJ'J 
3»* 
12  4 
-.2  1 
90  , 
92R 
9 A 9 

72> 

6 tn 
5 .3  S 

8 '(  O 

367 
2 93 
223 
3 9 | 
n5ti 
- 3 3 


4 a . > 3 2 

5 a • 6 6 a 
57*273 
♦ a . 3 3 9 
r:.2S9 
3 7 . j 7 f, 
3l.S73 
3 ' . 6 a 5 
23*292 

2a. 939 

2 3*  96^ 
2 . 8 3 ** 

I 9 . 22 l 
I a.  | 27 
l-.. 377 
I 3 • 3 9 j 

J 2*252 
l7.H39 

I 7 . 2** 
9 » 9 j A 
9.373 
9.4  12 
a • 5 q 0 
9.833 


1 5 . 36  I 3 
13.719  8 

11 .933  5 
1 J . 5 9 0 6 
9.373  a 
3.627  9 
7.772  1 
7.099  1 

®.527  1 
5.774  l 

5 . 9 rj6  2 
3.922  2 
8.239  3 
3.792  5 
3.390  7 
3«  1 2*  8 

2.9  1 2 1 

2.533  l 
2.314  1 
2.  m3  2 
1.S93  3 
1.581  1 

1.3S7  2 
1.397  3 

1.368  7 


0 

I 


s 
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SYMMETRIC  CONSTANT  Hll . 

heat 

t-t  L rt_a 

ruux 

» r 3 

PR. 

C»  .703,  OPP. DH/kC, TC»  2n0.C, 

r.iHMA*  1 

.800 

PROPPRTY  EXPONENTS.  VIS* 

.roc-  ,c  on» 

. r C c , r p * 

.oco 

— --  - ---  — , — — .. — - _ - - 

8X/DPE0 

TA/TI  TV. /TP  Nlinx  PRr/74 

A A / p A NUBA 

ERra/28 

l .pC-03 

38.46 

5.364- 

A8.73P 

1 0 

. 5 4 A 

1 ♦ 1 5-03 

37.1R 

5.0*3 

AC . *64 

9 

. 963 

l . 3 c- a a 

30.  AP 

8 . 7 AO 

S7.294 

9 

.207 

1 *SC-03 

?®.»c 

4 , *«r 6 

53. A 21 

A 

.650 

1.76-03 

26.94 

4 .01* 

A 

• G 7 5 

2 » PC-03 

25.  «2 

3.  Ama 

4 A . c 7 2 

7 

• 5 12 

2.30*03 

23. *5 

3 . A"  A 

4 8 » C A 5 

7 

.022 

2.  AG-03 

27.73 

3*841 

8 1 . A 7 3 

6 

.620 

3 « CC-0  3 

2W*Q 

3 . 2C  P 

39  .C5p 
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• 199  .Qfj 

1 » 00-0  l 

2 • Qo 

1 • 3*a 

8 * * 3 

1 . 0 1 3 

1.232 

1 A » 3 Q t 

1 .836 

8 

• 028*oo 

7*00-01 

3 • 30 

Ul75 

8*37 

1 .266 

1*131 

9*121 

! .553 

2 

• 102*01 

3 • OU-O 1 

4 • 00-0 1 

9 • 00 
5 * 00 

1 * 1 1 0 6 
1 • n7  2 

8*33 
8 • 3 a 

1.176 
1.13  3 

1 *C35 
1 *C60 

o.72g 

a . 5 2 q 

1*929 

1*364 

4 

7 

* 1 19*01 

• 063*0  1 

5*00-01 

6*30-01 

6 • 00 

7 * 00 

I • <3*2 
1 • n H 0 

8 • 2 o 
6 • 2 a 

1 . 1 I t 

1 . 1 0 1 

1*045 
1 *035 

4*399 

a • 3 2q 

1*327 
1 • 3q6 

1 

1 

• 1 23*02 

* 699*02 
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LAMINAR  f,  a 5 FLO  A BETAFFN  p A R A L I EL  PLATES 
S Y H M E TRjf  CONSTANT  A a l * MFaT  flux 


P P . 0 ■ . 667  . OPP.nH/KC.Tr.  7.0.  ramma«1.667 


PROPERTY  F XPONFnT*:’  V|S«  ,75O*r0N.  .750  . OP*  .000 


4X  /DpEO 

TP/T  ! 

Tti'/rB 

*URX 

FRP/74 

A A / R A 

nUPA 

F R B A / 2 4 

PO-P  t 

J 

• nc-03 

1 • CO 

. rt  7 

3 4 

.69 

5 • P c A 

1 • 0 2 9 

65 

. 1 3 1 

1 1 • q9B 

1 

. P23-C  1 

1 

• 30-03 

1-00 

• PAS 

30 

. a a 

5*120 

1*032 

57 

.6  39 

9 , 7 7 P 

• 0 9 5 - c 1 

1 

•7S-C1 

I • CO 

• 0 73 

27 

* I ? 

4 . 4 C 7 

t *037 

*0 

• 79n 

p.  pa? 

2 

.458-01 

2 

.30-03 

; .00 

. P P 7 

?4 

. I 1 

3,939 

1.091 

4 4 

.325 

7.946 

2 

.850-0 1 

3 

.00-03 

1 • 0 1 

. OR  7 

2 i 

.95 

3*  *20 

| .046 

39 

. ?87 

6.569 

3 

. 797-c  1 

3 

• ao-03 

1 .•  0 I 

. n P P 

l9 

.97 

3 . ? B 9 

1 *05G 

36 

. 1 8 a 

6.035 

3 

• <>Sc-0  1 

4 

• So-03 

t • 0 1 

• 1 0« 

j8 

• 79 

2 . ? 3 4 

1.0*9 

37 

,755 

5. 94p 

4 

■136-01 

5 

• * 0“  0 3 

1*01 

1 

* 1 1 A 

.6 

. a 4 

? Z 2 1 

1-0*9 

29 

.9  74 

9 • 960 

4 

• 6 3 « - q j 

A 

• AO-03 

1*01 

• 1 7 * 

1* 

’•?0 

2 .•  9 1 7 

1*0^3 

77 

.677 

4.565 

5 

• 1 * 3-C  1 

8 

. no-Oa 

1,^2 

• 1 33 

14 

. A 1 

7 • 3?  A 

1,067 

75 

• 47  1. 

4 • t P.  5 

5 

. 765-0  1 

l 

• no-07 

I • C 2 

• l **3 

13 

• 4P 

2 M 3 1 

1 *07  2 

23 

. 1 67 

3.786 

6 

.579-01 

1 

. 30-07 

1 .*  0 3 

. 

12 

• 34 

i • ? 3 ? 

1 .•  p 7 e 

20 

. 775 

3.37? 

7 

*714-0, 

1 

• 7 S -0  7 

i -o** 

• JAP 

T 1 

•79 

)•  7 4 j 

1 . 0 « 5 

1« 

.433 

2.966 

P 

• 289-0  1 

2 

• 30-02 

1 . 05 

• J 7 a 

1 0 

. 39 

I . c Q 7 

1 .•  0 9 1 

! * 

• c " 7 

2.644 

1 

• 1 C7*q0 

3 

.00-07 

1 , 06 

. 1 a 6 

9 

.73 

1 » 46P 

1.096 

1 s 

.029 

7.372 

1 

• 32  3*00 

3 

•ao-02 

1 *07 

. | P 9 

p 

• 35 

J . 3 p A 

1 .•  0 9 8 

1 9 

• P 9 3 

2-207 

1 

•501*00 

4 

•*0-07 

J *C9 

• l9  • 

P 

• 99 

1 • 3 7 2 

i • 1 0 0 

I 3 

• O8  1 

2-027 

1 

.762*00 

S 

• SO-07 

1*11 

• 1 PC 

a 

.75 

i • 2 7 n 

1*10  0 

17 

• 29n 

1,884 

2 

• 097*00 

6 

• AO-0? 

i • 1 3 

• 1 P 7 

a 

•AO 

1 . 2 7C 

1*0" 

I i 

I I 

• 663 

j.769 

2 

. 36  1 *00 

a 

.OC-O7 

1 , 1 6 

• ! e 2 

P 

• 49 

! , 7 P ? 

1.098 

. 09  1 

1.664 

7 

. 766*00 

l 

• no-o  1 

1 • 20 

•173 

A 

.43 

1 • 1 94 

1 . 0 9 9 

1 0 

.539 

1 . 56j| 

3 

.361*00 

2 

• no-o  1 

l • **0 

• 13  3 

a 

• 35 

1*1*7 

1 1 , Q 7 a 

9 

• 3 a p 

1-353 

6 

. 7 9 p *nn 

3 

. no-o  t 

1 • A 0 

• 1 n a 

a 

•33 

1,13* 

1 . r>A5 

a 

. 98  3 

1 .2  77 

1 

* o99*n  1 

4 

•00-01 

i • e a 

. p PA 

a 

• 3 1 

1 .•  1 1 3 

1-055 

a 

. 77P 

1-235 

1 

•607*0  1 

M M Ml 

"J  MM 

.3O 

1.048 

1 3 1 M 

.199*0  1 

5 

• P U • 0 * 

/ # C u 

.07? 

a 

1 • 09  A 

a 

• 657 

\ • * 1 0 

2 

6 

•no-o  1 

2 • 20 

•OA  1 

a 

• 79 

1 *099 

1*092 

8 

.565 

1.193 

2 

• 898*0  1 

\ 
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LMMIUmR  GMb  FLOW  bETwEEN  r AaALLc-L  PLATES 


5YMF<t.TRIL  COI(ST«TiT  'wALL  mEaT  1-LUa 

t 


Frt » u- 

• o o 7 * \jHH  • UH  / Ku 

• T 0—  1 uU  • • 

o h MN4  — 1 • 

60  7 

RHoFEH  T r 

t At  Of  jC-i  4 Tb  * Vi  b — 

. u o 0 * Com—  • 

OoO  »CP  = 

.000 

4A/LFLU  Tri/Ti  TV*/ I H IjUHA 

FRd/24  yA/bA 

NUH A FRHA/24 

PO-PI 

1 • 15-U3 

32.  b9 

5.  ibU 

o 1 . obo 

1 0 • lfj7 

1.860-01 

X • 3u-U3 

ju  . ab 

4 . a 2 u 

b7 • ooa 

9.516 

1 .979-01 

1 .bu-03 

2b. 97 

4.512 

b3 • 9ob 

8.861 

2.127-01 

1 • 75-u  3 

7 . lu 

4.21a 

bu  • 2o<) 

8 .224 

2.303-01 

2 • uu-u3 

2b . bo 

3 • 9bO 

47.209 

7.700 

2.464-01 

2*30-03 

24 . ua 

3.712 

44 . 341 

7.194 

2.649-01 

2 • ou-u  3 

22. d4 

3 • 49a 

41 .ob2 

6. 780 

2.820-01 

3 • ui/*uj 

21.31 

3 . 289 

39.298 

6.327 

3.037-01 

3 • 3u~03 

2l> . oo 

3.157 

37.642 

6 . 044 

3.191-01 

i.ou-uj 

19.92 

3 . u3b 

3o • 19b 

5.798 

3.34U-01 

4 .uo-oi 

19. ua 

2.066 

34 .52b 

5.514 

3.629-01 

4 . 5U~U b 

la. 19 

2.731 

32. 7ol 

5.212 

3.753-01 

| 5.uu-u3 

17.44 

2.O09 

31 . 2on 

4.958 

3.9b6-01 

5.5u-u3 

1 o • 7 9 

2 • b02 

29.979 

4 . 739 

4.170-01 

o . uu-u3 

lo  .2b 

2.40b 

20 .Hb7 

4 .54H 

4.3b7-01 

, O.tJU-tJb 

lb . 04 

2 • 3 0 4 

27 .602 

4.348 

4.592-01 

5 7 • 2 0 - u 3 

lb. 13 

2.224 

2o . 657 

4.174 

4.809-01 

6 • 00“ Ob 

14  . b4 

2.119 

2b . 474 

3.974 

5.087-01 

9.0u-ub 

13.93 

2.014 

24.22b 

3 • 7 fy  1 

b. 415-01 

r 

1.0U-U2 

13.41 

1.931 

23. 1 70 

3.582 

5.731-01 

1 • 15-U2 

12.78 

l .aia 

21  • 8bo 

3.358 

6.179-01 

1 • 30-U2 

12 . 2o 

1.732 

20.779 

3.175 

6.604-01 

■ 1.50-02 

11.70 

1 • o3o 

19.6U5 

2.975 

7.141-01 

* 1.75-02 

11.1b 

1.543 

10.437 

2.77b 

7.774-01 

2 . 00-02 

10.72 

1.469 

17.499 

2 • b 1 7 

8.375-01 

- 2 • 3u  — 02 

10.31 

1.39b 

lo . 508 

2.462 

9.0ol-01 

4 2*60-02 

9.9o 

1 .34U 

lb. 844 

2 .336 

9.716-01 

3 . Uu-02 

9 . 04 

1.279 

lb. 039 

2.148 

1 .055  + 00 

3 • bU“U£ 

9.43 

1 .242 

14.538 

2.113 

l . llb  + 00 

, 3 • bo-02 

9.2b 

1.209 

14 . lob 

2.0  34 

1 .175+00 

4.UJ-U2 

9.U7 

1.172 

13.611 

1 .954 

1 .251+00 

4O0-U2 

a.a9 

1.13b 

13.  [)9b 

1 . o65 

1 .343  + 00 

m 5.UU-02 

a . 7b 

1.10  7 

12 • ho8 

1 .74() 

1 .432  + 00 

► 

| b.5u-U2 

<1.04 

1 .Ufib 

12.307 

1.727 

1.520+00 

b. UU“02 

a.bb 

1.060 

11.997 

1.673 

1 .b0o  + 00 

o « ou-02 

a. 4a 

1 .051 

11 .hhl 

1 . u 17 

1.708+00 

7 • 20-u2 

a. 42 

1.039 

11.411 

1.564 

1 .aua+on 

d • uu-u<: 

a . 3o 

1.027 

11.109 

1.51b 

1 .940+00 

9 . uo“U2 

a . 32 

1 • u 1 7 

lo .402 

1.461 

2 .103+00 

1 .Uu-ui 

a. 29 

1.012 

10.552 

1.416 

2 .2ob  + 00 

1 • 50-u 1 

b. 2b 

1.003 

9.741 

1 .274 

3.U69+0O 

2.00-Ul 

a .24 

1 .U01 

9.404 

1.204 

3.870+00 

2 • 30  — 0 1 

6.24 

1.001 

9.1/1 

1 . 168 

4 ,b7l+0o 

3 • Ou-O 1 

a .24 

1 .001 

9.01b 

1 . 140 

3.471+00 

u 

3OU-01 

a .24 

l.uoi 

a • °o4 

l .120 

n.272+00 

4 *00-0 1 

0.24 

1 .001 

a • 8 2 0 

1.108 

7 . 072+00 

4. 5u-ul 

d .24 

1.001 

a • 755 

1 .093 

7 .M73  + 00 

b • Ou-u 1 

0.24 

1.001 

0.7o4 

1.044 

a .0 73  + 00 

5.30-ui 

6.24 

1.001 

a . ho  1 

1.077 

9.4/4+00 

b-ou-ul 

a .24 

1.001 

a.626 

1 . 07U 

1.027+01 

o • 5u-u l 

a .24 

l.uOl 

a • 696 

1 .065 

1 . 107+01 

- * 
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L A*  In  Aft  p*  s PLO'V  BET.tCEN  P»9alLEi.  pLATES 
SYlMETRfC  ConSTjNT  * *1-1.  Mp  A T Fl_.»X 


PPiO»  .467.  3PP . OH/ <0.  T*»»  -2.0,  GA-*A»J»667 

p«opepty  exponents  , vi<;»  . 750  • con»  .75 A , cp»  »ooo 
tX/OPEO  T0/t1  T * / T 9 NUB*  PPr/ ?4  *4  / qA  nUBA  P*a a7  24  PO-PI 


0 0-03 
1&-03 
30-03 
50-03 
75-03 
0 0 - 0 3 
30-03 
40-03 
00-  0 3 
30-Q3 
60-03 

0 0-  0 3 
50-03 
00-03 
50-03 
00-03 
60-03 
20-03 
00-03 
00-03 
00-02 

1 5-Q2 
3O  -0  2 
50-02 
75-02 
00-02 
30-02 
60-02 
00-02 
30-02 
60-02 
00-02 
50-02 
00-02 
50-02 
00-02 
60-02 
20-02 
00-02 
00-02 


3 3 

.9H2 

34  . * 9 

5 *2  62 

.971 

00 

.938 

32*57 

4.898 

.9*9 

03 

• 9 35 

30*  05 

4.634 

.967 

00 

• 93  1 

28 . 97 

4.267 

.9*5 

Co 

• 9 2 6 

27*09 

3.955 

.9*3 

03 

.921 

25.54 

3.733 

.9*1 

00 

.916 

24.0c 

3.455 

. 9c8 

99 

•912 

2 2-97 

3-252 

.956 

99 

.906 

2 | . 47 

3.057 

.953 

99 

.902 

20’  6 7 

2.936 

.95  1 

99 

« a 9 8 

1 9 . 8 a 

2.766 

.949 

99 

• 8 9 3 

1 9 . 0 3 

2.634 

.947 

99 

. a 8 3 

1 8 . 1 4 

2.494 

.944 

99 

.083 

1 7 . 3a 

2-366 

.942 

99 

. a 7 8 

1 6 • 7 -j 

2.258 

.939 

99 

. a 7 4 

| 6*  1 7 

2*165 

.937 

99 

. a 69 

1 5 » 57 

2-064 

.935 

9? 

.964 

1 5 • 0 * 

1.984 

.932 

9fl 

. a 5 8 

1 4. 4? 

1.878 

.929 

93 

. 0 5 | 

1 3 • 9c 

1.775 

• 926 

9a 

. a 4 5 

1 3*  3i 

1 .687 

■ ’23 

9a 

• a36 

,2*69 

1 .579 

.9  1 9 

97 

. A28 

1 2 • 1 7 

1 .494 

.9,5 

97 

.818 

1 1 • 60 

1 .394 

• 9 i 0 

9 7 

• 807 

1 1 ♦ n c 

1 . inn 

. 9 -.5 

96 

.797 

1 1 T 
l C • 6 n 

l .226 

• 9 A 1 

95 

.787 

1 0»  1 9 

1.153 

.896 

95 

.777 

9.8c 

1 *09  1 

.89  1 

94 

• 765 

9.5  n 

1*027 

• 886 

93 

.757 

9.2a 

.989 

.883 

93 

• 7 *♦  9 

9 . 1 1 

.953 

.879 

’2 

.743 

8.9/-. 

.9(19 

.875 

9 1 

.729 

8.7) 

.86® 

.87  1 

’0 

• 719 

8 • Sc 

.8  35 

.867 

8 9 

aa 

.7  0? 

•7  00 

8.4a 
8 . 3 1 

.8o7 

.792 

.8*3 
.8  59 

8 7 
86 

• *89 

• *73 

8.24 

8>u 

.757 

.735 

.955 

*852 

84 

• *64 

ft. 90 

.7o9 

.847 

82 

• *47 

a . 0 1 

.68  1 

.841 

65 
6 1 
57 
5 7 
5; 
*7 
*4  u 
8l 
3 o 
37 
3a 
3 a 

3? 

3i 

29 

20 
27 
2a 
2c 
24 

21 
2 7 
2.; 
1 9 
t 0 

1 7 

t A 
1 5 
15 

1 U 
14 

13 
1 7 

l? 

1 7 

l I 

l I 
1 » 

1 7 

I r> 


1 39 
075 
6 9 0 
989 
234 
29  1 
362 
*S» 

3 14 
658 
211 
54C 
773 
277 
989 

8 65 
689 
663 
48o 

2 3 0 

174 

858 

780 

605 

336 

497 
535 
841 
035 
533 
099 
603 
q37 
6 S 6 
292 

9 8 0 
660 
387 
0 79 
764 


I 0 

9 

9 

A 

7 

7 

6 

6 

6 

5 

5 

5 

4 

4 

4 

4 

4 

3 

3 

3 

3 

3 

2 

2 

2 

2 

2 

2 


• 5 4 l 

1 

•647-01 

.818 

1 

•761-01 

.24  1 

1 

•870-01 

.598 

2 

• 002-0  1 

.956 

2 

•155-01 

.440 

2 

•298-01 

.937 

2 

•856-01 

.524 

2 

• 603-0 1 

.07  4 

2 

♦787-01 

.796 

2 

•9 17-01 

.549 

3 

•039-01 

.266 

3 

•193-01 

.966 

3 

•375-01 

•715 

3 

•547-01 

.49  8 

3 

• 7o8-CI 

. 3 Q9 

3 

•962-01 

>111 

4 

•036-ot 

.939 

4 

•2o2-o  1 

.740 

4 

• 4 1 1-01 

.S30 

4 

• 653*0 1 

.353 

4 

• 982-01 

•131 

5 

•197-0, 

.95, 

5 

.49,-0, 

.754 

5 

• 849-0  1 

.557 

6 

• 257-0  1 

.400 

6 

• 628-01 

.247 

7 

•036-0  1 

• 123 

7 

•409-01 

.987 

7 

•860-01 

• 9 0 3 

a 

•173-01 

.830 

8 

•465-01 

.746 

8 

•822-01 

.657 

9 

•230*01 

.593 

9 

• 602-0  1 

.520 

9 

•942-01 

.466 

1 

• 026*00 

• 410 

l 

•060*00 

. 362 

1 

• 09 1 *00 

.307 

I 

♦ 1 28*00 

• 250 

1 

• 169*00 
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LA-I*aR  G 4 S ftO»  9 E T ' £ E N P 4 9 ALl  EL  PLATES 
SVMMgT^if  CONSTA'IT  •»  a l L jT  r |_  U X 


?S,Q.  . 9 q 3 , OPP.TW/<o.TC»  |n3.3.  6AMM*»|.**7 


PROPERTY  fXPOMfwT*,'  YfS*  .75n'foN»  . 750  , O'**  .000 


9 l /SPED 

T9/T1 

T 6/  T9 

»:)3  t 

599/74 

4 4 / 5 A 

N'J9  4 

P98A/24 

PO-PI 

l .O0“03 

1*10 

3.296 

37.67 

19.071 

2.  1 77 

25.169 

4 . 635-C 1 

t • 10-03 
I . 75-0.3 
2 . 30-03 
3 . n o-0  3 

1.13 
1.17 
1 .23 
1 . 3 0 
1.35 

3 . ^9  A 
3 * t;  q ? 

3 3 . * 9 
?9  . *0 
26.36 
7 3.95 
2 1.66 

1 1.993 

10.313 

9.340 

9.193 

7.3*0 

? * 27  1 
2.379 
2.459 
2.574 
2.555 

6 1 . o9o 
5 3.  1 9 5 
46.740 
9 1*237 
37.9U 

22  . 299 
19.050 
[6.562 
19.395 
13.055 

5.7  15-01 
7.184-01 
9.997-01 
1 • 1 2 6*00 
l . 323*00 

3 • A *4  A 

3 * A Q 7 

3.  *0-03 

3 . 7n0 

9.50-03 

1 .95 

3.6*5 

1 9.  *2 

6.779 

2.576 

33.97? 

11.599 

1.61 2*00 

5.50-03 
* . A 0 - 0 3 
8.O0-03 

1 .55 
1.66 
1.90 

3.59* 
3 » 5 ft  t 
3.3*7 

17.91 
l'  6 . 9 9 
1 5 . | 1 

6.755 

5.77n 

5.445 

7.577 

2.560 

2.522 

30*915 

24.(71 

25.593 

1 0. 3*7 
9.442 
9.5  33 

1 .937*00 

2»3o9*qo 
2 • 78  1 *00 

1 .00-02 

2. 00 

3.179 

13.67 

4.999 

2.452 

22.973 

7.625 

3.436*00 

1*30-0? 

2.30 

2 • 9 n 7 

12.23 

4.5*9 

2.329 

20*006 

6.714 

4.61  0-00 

1.75-0? 

2-75 

?.5*9 

10*90 

9.0?0 

2.196 

17.194 

5.837 

6.447*00 

2. 30-0? 

3 .30 

2.739 

10*00 

3 . 4 ?4 

1.95: 

15.015 

5.109 

8.913*00 

3 • oo-Q  ? 

9.00 

t .9«  t 

9.39 

2 . « 1 1 

1.753 

13.272 

9.4U 

1 .232*0 1 

3.50-0? 

*4 . 60 

1.759 

9.12 

7.4|* 

1.623 

1 2*  2*4 

3 • 9 3 1 

1 . 545*0 1 

9.50-0? 

5.50 

1.573 

9.99 

2.020 

1-942 

1 1 .290 

3.373 

2.099*01 

5.53-0? 
*. 50-0? 
* . n o-n? 

J • (30-0  i 
2*00-0  1 

6 .50 
7 .60 
9.00 
1 1.00 
7 1*00 

1.932 
1 . 3 .3  3 

1.751 
1.179 
1 . 0 5 4 

9.74 

•» . * 3 
9.53 

9.45 

9*31 

1.734 

1 . S 9 7 
? . 4 r 5 

1 . 3 4f* 
1.214 

1.375 
1.294 
1.226 
. . 7 * •» 

10.575 

2.919 
2.559 
2*296 
1 .96* 
1.599 

2 . 66  2 *0  1 
3.9lo*ot 
4.494*01 
6.369*01 
2.474*02 

10.099 

9.594 

O . t BC 

1 • 

* . <4  l «4 

3*00-01 

3 1*03 

1*030 

4.29 

1*137 

1 *a?9 

8.2Q5 

1 . 46o 

6. 194*02 

9* no- at 

9 0 • 9 9 

1 • 0 1 9 

9 . 27 

1*113 

1 ■ n 1 3 

4.119 

1.394 

1.273*03 

5.00-0  1 

5 3.99 

1 • 1 1 ? 

4.76 

1 * M 9 

1.012 

4.063 

1 • 3 7 0 

2.363*03 

5*00-0  1 

6 3 • 9 8 

1 . an  9 

4.79 

1*297 

1 *009 

4 • ri  3 3 

1.367 

9.29  , *03 

t=A*  I. 4*9  «>*>  fLO“  ftE-T.-rETM  ^ ALL  Ff  rLAT€S 


> y 1 4 £ T ^ j ^ C 0 *4  S • * •*  T rV  A L L M p A T r i,^.  A 


PD 

t r.  x ,4f 

C«  wF'P.O*'*-'?. 

Tr.rn 

1 0 

t 

• 6 A 7 

pSnPERrY  £ x Po^r^T  5 , v|<;« 

. 7S~  t 

cOm*  • 

75-  , CP* 

♦ 0 CO 

u * / D p E 0 

ra/r  j 

Tn/T5  ;;uR*  F B 3 / ? 4 

Vi 

i/pi 

**  y B A 

Fsp  A/24 

PO-P  1 

1 • C3-03 

1.0! 

1 . ? 6 7 J A . ? 4 

n . 6 - 7 

1 

.13“ 

6 a . 3 0 l 

15.603 

I 

7 15- j! 

i • 3 0 - u 3 

» • 0 l 

1 • 2 9 p i 2 • 7 ■* 

7.60“ 

1 

•15. 

6 » • 2 ^ ^ 

13*975 

2 

0l6-ul 

t .75-33 

1 * v 2 

1 • 1 3 7 J e • 7 a 

6 . s :3  4 

1 

• 17:.. 

51.37? 

12.135 

2 

*?5-'l 

?. 33-33 

i • I 2 

1*37,  25 . 5S 

A . / 44  fl 

1 

.1°: 

4i  . j5  1 

1 -. 744 

2 

087-01 

i • o ; - o 3 

1 • .3 

1 • u j fe  ? 2 • 9 * 

5 . 4 a 2 

3 

• 211 

ST  .46-1 

0.545 

3 

*29-0  1 

3 • * 0-33 

; .ci« 

1.44*14  2 1 • 1 A 

5 • - 3 9 

1 

.22ft 

3 0 • 3 “ 3 

« . 9 1 3 

3 

57  1-3  1 

u • 5 0 - 0 3 
5.5;-o3 

1 * 0 4 

I>v4 

l • *4  7 *>  t V • 3 3 

1 • =4 1 :*-  1 7 * 8 a 

4.420 
4 • 2 4 7 

i 

i 

• 245 

• 2 4 2 

3^.667 

37*609 

7.999 
7 • 3 l 7 

4 

5 

4 9 8-01 

164-31 

3.947 

J 

.273 

.294 

20*226 
2a  * 657 

6*753 

5 

P63-01 

721-01 

tt.CO-Oi 

i «ca 

1 it66  i S » 4 6. 

3.663 

1 

6 . 20 7 

6 

t • •?  3 “ 0 2 
i » 3j  - 0 2 

1 • 1 o 
1 • 1 3 

| • 5 9 4 j *4 . 2 u 
1 • A2  ! 1 3 • Cr- 

3.356 
3 - 4,i  4 3 

1 

1 

.31! 
. 3 3.j 

2u  . 37? 
2 » * 0 0 0 

5.629 
5. *>21 

7 

9 

9o5-at 
A 2 l " 5 1 

, .7&-02 
7 • 3 0 - 0 2 

1 ♦ 1 7 
l • 23 

| ♦ A 3 0 1 I ♦ 9 1 

1 * — ‘*•1-  i 0 • - “ 

2.743 

2.4*7 

1 

1 

• 3 •«  5 
. 3 S 3 

10*262 

17.247 

4.4  i5 
3 • 9 2 7 

1 

l 

2 1 2*00 
5 1 •*  ♦ 0 n 

1 *.~0-02 

1 • 3 -i 

1 • a 2 3 1 ; » 1 “ 

2.265 

1 

. 352 

I * . 5 4/, 

3 • 5r  9 

i 

3 97-Ct-, 

3 • 6 > - *)  2 

1 • 3 o 

l • AOO  9.7, 

2.150 

1 

.346 

1 4 * 5 1 7 

3 * 2 5 0 

2 

2 2 8 * 3 0 

« • 5 ,-92 
5.  S0-.32 

l »“5 
l *55 

1 ♦ * » f V • 3 1 

1**13  9.2a 

7 • yf.3 
1 . a75 

1 

1 

.33  1 
* 3 1 2 

1 ■»  • ‘'CC 
17*526 

2*962 

2*727 

2 

3 

735*30 
3 t 7*co 

a * 6 O - O 2 

I .4* 

1 # a 3 a « a o 

1.770 

1 

* 23  9 

1 7 * - 3 3 

2*532 

3 

9fl  1 *00 
0 4 9 *00 

0 • o>-o2 

l * 4. 8 9 . 7 A . 

1 • 6 6 7 

1 

• 242 

1 7.*  2u 3 

2*344 

4 

» ♦ <3  4.-0 1 

,.cc-o« 

2*0:'. 
3-r : 

1*343  a . 6 a 
1 . 1 73  I . 4„ 

5 . 5 A 9 
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1 

1 
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• 1 3^ 

1 , • S i 9 
c • 2*,C 

2.  1 47 
1 .721 

6 

l 

22°*on 
6 2 9 ♦ 0 1 
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«4»0w“0l 

**  • Cv 

- • £.!V 

1 ♦ 1 06  8 ♦ 3a 

l'r72  0 . 3 1 

1 . 2 | 5 
i • l 5 7 

1 

1 

• 035 
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4 ♦ 3 3 7 
Q • 6 Q 7 
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1*93» 

2 

4 

8 5 3 * 0 1 
7 1 7 ♦ 0 1 

*♦00-01 
a* 90-3 1 

*>  • 00 

7 • -o 

i^r-.sz  e . 3 , 
1 * . 1 5 0 3 * 3.’. 
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i 

1 

.•0“5 
• 03  5 

a • 

0 • 3 7 5 

1*302 
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7 

1 
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LAMINAR  f,  a S F I 0 Ji  H F T 6 F F N PARALiFL  PLATF5 
SYMMETRIC  CONST  ANT  .'/ALL  HEaT  FLUX 
P R . 0 « • 2 C 3 i OPP.DH/KC.TO*  103.0.  c,  AMMA*1.667 
PROPERTY  EXPONENTS.  V i S ■ . 7SO^CO‘i*  .750. CP*  .033 


SX/DPEO  TB/Tl  T.V/TO  NUBX  Fr  B/7S  .*!  A / 9 A NUBA  F R R A / 2 “ K 

1. 00- 01  1 • 1 G .1.170  SO.  10  1*7. 5i7  2.  1C9  7s.  1 57  3S.261  .25B 

1.10-01  1 • 1 3 3.765  35. bo  I6.PSA  2.7nO  A 5 • <4  7 i 3o. 20*  .310 

1.75-03  1 . 1 7 3 . **  1 0 11. '4  5 13.879  2.303  5 6 . a 7 p 7 5 . 9 b 2 . 3 B I 

2.30- 03  1.23  3.5rts  27.91  11.7«o  2.380  S9.9SS  22.39s  .S60 

3.00- 0  1 1*  30  3 . 5 a 2 7 4.76  9.697  2 .SS7  S3. 987  t 9 . [ p s ,5SR 

3.60- 03  1.3*  3 . 5 a 8 27.80  8.716  2.S79  so.  773  17  . 1 32  . 6 1 9 

5. 50- 03  1 • u5  3 . 5 s 0 2 3 *59  7.9 26  7.503  3A.J15  | S . 9 ( s .716 

5.5  0-03  1.55  3 « S 7 5 1 8 . 7 8 7. 3 7 7 2.50**  32.700  13*222  .817 

4.60- 01  1*66  3.39s  1 7.7s  6.688  ?.S»3  29.85s  11.83)  .970 

8.00- 03  l.so  3.?a5  15.79  6.503  2.S56  27.069  I0.*05  1 . c S5 

l.nC-O?  7.00  3 .nas  is.  76  6. Cue  2.390  2s.|39  9 .h05  1.712 

1.30- 02  2.30  2.817  17.72  5.7t6  2.276  71.071  8.279  1.SS5 

1 . 75-02  2 • 75  2 » 5ps  n.33  5 . I a ( ? . 1 q 3 I8.061  7.27s  1.771 

2.3Q-07  3.30  2.192  (0.39  «.5r>6  |.9i5  15.736  6.S60  2.128 

3.00- 07  s . OC  l • 9 p 5 9.77  3 « 7 n 5 1.72S  1 3.886  5.655  2.51  Q 

3.60- 07  a.60  1.778  9,cO  3.IAI  I.598  1 2 . 8 s 1 5 *068  2.77s 

S.5C-0?  5.50  l • 5 s 7 9.77  2.610  1 . s 6 2 J 1 . 7 9 S S.36Q  3.Q87 

5.50- 07  6*50  l.S(*  ».n8  2.203  1.361  ll.n38  3.7as  3.35i 

6.60- 07  7.60  1072  8.92  1 . 9 1 5 ) . 2«  5 (0*s73  3 * 279  3.57Q 

8.00- 07  9.00  1.7SS  8.76  1.687  1.720  9.975  2*836  3»78i 

1 • no-0  1 11.00  1.(75  b _.  a J 1.5?A  1.(60  9.505  2 • '4  | 8 s ,po3 

2.OQ-0I  21.00  1.058  8.36  1.766  l.jjSS  8.581  |.76S  5.2S8 

3.00- 01  3 | • 00  1*030  ".11  1.158  J.02P  8.308  [.5S5  6.J75 

8.00- 01  Sfl.99  1.CI8  8.78  |.M6  l . Q 1 8 8.186  1 . S S 7 7.367 

5.00- 01  Sp.99  1.01?  8.75  1.097  1.012  «.|16  t.397  9.286 

6 • o 0 — 0 1 *0 • ^ ^ I • Cp9  8.78  1.q9S  I » go9  8 . p 7 2 1*  368  |2.*c7 
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laminar  gas  eloa  between  Parallel  plates 

SYMMETRIC  CONSTANT  *ALL  HEAT  ELUX 
PR.O»  *200.  OPP.UH/kO.TO»  SO. Or  GAMMAs | .667 
PROPERTY  EXPONENTS,  Vis**  .750*C0n«  .7SO,C Pm  .GOC 


HX/0PE0  Ta/TI  TK/TB  N U 8 X ERB/21  «a/8A  NUBa  ERBA/2h  K 

1*00-03  l»OS  2 • 1 S 2 39.8s  17.196  l.59o  71.627  2 9.129  *188 
2*00-03  1*  10  2*137  29*15  12.322  1*753  53.831  21  *622  *301 
3*60-03  1*18  2*636  22.87  8*981  1.886  10*889  16*210  *116 
6*00-03  1*3o  2 • 7 1 5 |8.l2  7*002  1.969  32*  182  12*138  *627 
1*00-02  l*5o  2*652  11*89  5.671  l*99i  25*283  9.177  .88! 
2*00-02  2*00  2*301  ll**<2  1*5q8  1.867  18.191  6.7q5  1*1o7 
3*60-02  2*  80  1*811  9*77  3.396  1.622  13*96o  S«o72  2*076 
6*00-02  1*00  1*182  9*17  2*381  1.385  11*516  3»8ol  2*769 
1*00-01  8*oo  1*218  8*77  1*711  1.212  lO*05o  2*731  3*152 
3*50-01  |8*5q  l*n36  8.35  l*2l5  1.q31  B.lli  l*6o5  7.tn7 


L AM  1 NA R GAS  PlOR  8ETAEE*  PARALLEL  PLATES 


SYMMETRIC  CONSTANT  RALL  MEAT  ELliX 
PR  *0*  *200 , QPP.OH/kO.TO*  25.0*  6AMMA»1*667 
PROPERTY  EXPONENTS,  VlS»  *75O,C0n»  .750. CP»  *000 


IX/DpEO  TB/Tl  TR/TB  n'uBX  pR»/2l  Ra/BA  NUBA  FRBA/21  K 

1*00-03  1*03  1*602  39. 8n  11.757  1.3qS  75*001  2S.273  *112 
2*00-03  | • OS  |.78l  29*29  i0»*87  l.lot  51*111  18.782  *223 
3*60-03  1*09  1.915  22*75  8.35’  1.193  11*212  11*156  .326 
6*  00“03  1*  15  2*o62  l8.1i  6*57o  1 *568  32*59i  ll*1o7  *156 
1*00-02  1*25  2*121  15*09  5*210  1*623  25*813  8.9i1  *637 
2*00-02  1 • 5o  2*0 **9  ll*77  1*01*  1 *627  18. 977  6*359  |*008 
3*60-02  l*’0  1 *813  10*00  3 * 2 1 3 1 *532  11*787  I.81I  1*188 
6*00-02  2 • So  |»S13  9.27  2*508  1*388  12*231  3*759  2*0*2 
1*00-01  3*5q  1*311  8.89  l*9o®  1*211  lo«565  2*867  2*7o9 
3*50-01  8*75  | » q56  8*37  l*2l>  1*q5o  8*6|2  1*61Q  1.868 


Inlet  pressure  was  relatively  low  on  these  two  calculations,  leading  to  Mach 
number  greater  than  0 .2  for  x*  > 0-3. 
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t.  AM  I N AH  SAS  FLO#  BETWEEN  ^aRaLLEl  PLATES 
SYMhETRjC  CONSTjMT  maLL  MEAT  ELu* 

pR»0»  *200,  «PP . OH/<0. Tn«  lO.O,  SAmMA»1.667 

property  Exponents,  v I s»  . 7S0 . cOn«  . 7 5<7 , cp«  . ooc 


| mx/DpEO 

* r—  - ir’-_ 

T8/TJ 

t«;ts 

MUR  X 

FRa/2* 

A*/RA 

luliaA  ERqa/24 

PO-P  l 

* . 00-03 

1 • 0 t 

t . 7*6 

3*. 

0* 

12.123 

3.12** 

7^  .3*6 

2 1 .958 

1 • 266-0 1 

> * ji 
» *01 

1 . ?7i 

i I » 30-03 

35. 

**i 

» 3. *79 

1*139 

6 a . 5 * 1 

19.518 

1 .*95-01 

li  7.7&-03 

l *02 

t • M 3 

31  • 

07 

9.6  n7 

1.158 

57  .*35 

17.107 

1 • 8 1 3-0  1 

?. 30-03 

1 .02 

J * 1&  3 

27. 

*9 

8.633 

1.177 

sT.ao* 

15.152 

2.171-01 

a 3*00-03 

1 *03 

1.388 

2 * • 

5 r» 

7.629 

1.197 

*n. 105 

13.965 

2*596-01 

5 3 • 60-03 

1 O* 

1 • « l 5 

22. 

6 * 

7.083 

1.211 

*7 . *7o 

12.*  1* 

2* ’*0-01 

*.50-03 

1 *05 

J *9*9 

20* 

6 7 

6 . 3 6'  8 

1.22* 

i7.**i 

1 1 .2*3 

3**33-01 

5.50-03 

1 *04 

1 *47* 

1 a • 

*3 

5.89S 

1*2*6 

3*  . 1 7 6 

l 3 • 2 7 6 

3*’59-ci 

4.60-Q3 

1.37 

t * 507 

t 7. 

6 c. 

5 . *03 

1.261 

3 7 * * 7 7 

9. *6* 

* • 5o3-0 1 

*»00”03 

1*08 

l *OH 

(6. 

37 

4 . 9 9 0 

1.277 

2a . 8 8 0 

a. 676 

5.  18  1-01 

» *00-0? 

1 • to 

| . P.62 

|S. 

04 

4 . 5 1 * 

1 . 299 

24. 185 

7.835 

6*111-01 

a 1/30-02 

t • 13 

1*58’ 

1 3 * 

7 n 

H .058 

l • 3 1 3 

2a»  3*c 

6.99* 

7**56-01 

| l « 7S-g2 

t . 17 

1*438 

1 2 • 

*7 

3.6o2 

1.32* 

2A.56J 

6.056 

9.9,5-0, 

7.30-02 

l *23 

» *6ll 

1 1 • 

39 

3.2** 

1.337 

la • 355 

5.3*0 

1 • 1 76  *00 

3*  oO”0 2 

t * 3o 

1 .5*7 

1 3* 

3a 

2.937 

1.337 

1 A • * * 6 

*.727 

1 * *79*go 

J 3 • 6 0 - Q 7 

1 • 36 

1 • 57  6 

JO* 

i 3 

2.753 

1 • 332 

15*369 

9.3*5 

l . 72**03 

* • 50"02 
5*50-0 2 

» *HS 
1*55 

J • 4 38 
1 » «9  4 

9. 
9 . 

69 

H <5 

2.596 
2*37  l 

1 .31* 
1.3-30 

I*. 1*8 

l a*  I9  l 

3 • * 30 
3.59  , 

2 • 1 | 8 ♦ 00 
2 • S 6 3*  00 

a 4*40-02 

1*66 

1 * * 1 7 

9 . 

2 7 

2*223 

l • 27* 

17. *39 

3«3  10 

3»  060*00 

| a.qO-02 

1 • flG 

1 * jJS 

9 . 

34 

2.37  7 

1.25* 

l 7 . 7 *♦  2 

3*0*0 

3*7  2 1*00 

1 *00-01 

2 • on 

l * 333 

8. 

*3 

l .925 

1*222 

1 7»q46 

2*757 

9 » 7 , 6*03 
1 .085*01 

7*0C-0J 

3 . 00 

1 M 7 0 

9 . 

l .5  l * 

1.128 

9.638 

2*098 

3«oo-oi 
* • 03-01 

*•30 

5 • on 

t • 1 3* 
l *.-»7  1 

a * 
9 • 

4 4 

*A 

t • 3 l 9 
J .225 

I.38* 
1 • -5  * 

«»07  l 
a.  797 

| .782 
1 .615 

1 • *03*0 1 

2 • 9 5 3 ♦ 0 1 

5*00-01 

A. 00-01 

6.30 
7 • 00 

1.-.52 

8 . 

34 

1.172 

1.0*5 

9.612 

1 - S 1 5 

*.3«9*o. 

1 • n*0 

a . 

33 

l . 1 * a 

1035 

a . 9*3 

1 **50 

6 • 168*3, 

kk  - 


I 

I 

I 

I 

I 

I 

I 

I 

I 


lA/-I‘.aR  r,  AS  (l5*-  4 F T a £ F f*  PaRaLLEi  pLaTfs 
SyMm^TKiC  CO^ST.wT  -‘•ALL  HE^T  fU-i 


P»tis  . 2 • LPP.Dw/vO.Tr,*  2 . ~ , GAukA*  I «Af7 

pSopr Rr>  f<Pr.Xp.‘<T5,  v|ci  ,7Sf  »fCt  * .75/--tfP=  •COC 
uX/pp^O  tB/Tj  t*,Ts  ajUBX  (-a/b a uUbA  fppa/24  Pc-PJ 


i .0*7-0  3 1.00  l.rtSc-  “O-.O-  7 1 ;1  . **0  l 1.02  5 7^,597  1 $ . 9 1 2 9 . 996-02 

It3u-C3  l«0c  l*.»56  35.57  0.1*47  1.026  6 i, . P 3 2 I 7 . S ^ S 1 • I — C.  1 

7.?5-c3  »*CD  1*'6h  3i.fi  P.:_65  j . 3 7 59. 165  15. 25*  1*35^-01 
7.3C-Q3  1 . JC  l • A 7 2_  2 .7  • S ».  7 . Ij  6 C 1 - C 3 <!?  5 7 • 2 2 n 13*397  1 . 5 9 r - n 1 

TC'J-03  »*Cl  !•*>»{  A.2R5  i - c-  *♦  1 H^l-j  1 l - ^ 7 C 1 . &3S-G  X 

3.45-03  I.'I  1*307  7 2 • 6 - S.777  i.f»H  t,  .669  ir.«37  2 »C32**Q'l 

o.5l-03  l*5i  l.*>9<i  2u»57  S . 2 1 **  l.C-VR  37.633  °.76f>  2.3JI-H 

5. c^-03  l * C i 1 • i C“  I4*9?  *1  7 3 4 i»05  2 3**.365  R.o«P  2 *694-0  1 

A.6C-03  l * C 1 lull  ,7.57  4 . 3°9  1 » C-  5 6 3,  • 663  6*156  2 ■ PR  3“C 1 

a.CC-03  1.02  l*il9  1 6 . 2 a «.r.c3  I.C-6C.  29. o67  7.  “56  3.23u"Cl 
7»nO-o2  l » C 2 1 ♦ i 29  i H . 9 a 3 .6“5  i . r 65  2a*3&6  6*7i5  3 .4A6-q  1 

T • 3 G " 3 2 l * C 3 l'l'l  *3*^  3*23°  1 • 0 7 l 23**4l  5 • 9 H 1 ‘,,32rci 

7.75-0  2 1 *03  | » 1 5 3 7 2 * 3 ^ 2.661  l‘C-7b  27.762  5.17*4  5*l?i"01 

? * 3 l-  ” Q 2 l*o5  1M64  ii*3  ..  2*-5?  1 » r.  8 H ic.ecc  *4  • S fc  4 6«i77-.2i 

3*0t“02  I • C 6 t * j 7 2 mj.50  ?.3ftR  | - c “ 9 I**767  R*045  7,35i~oi 
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